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Abstract 

Based on our recent NNLO dynamical parton distributions as obtained in the 'fixed 
flavor number scheme', we generate radiatively parton distributions in the 'variable 
flavor number scheme' where also the heavy quark flavors (c, 6, t) become massless 
partons within the nucleon. Only within this latter factorization scheme NNLO cal- 
culations arc feasible at present, since the required partonic subprocesses are only 
available in the approximation of massless initial-state partons. The NNLO predic- 
tions for gauge boson production are typically larger (by more than la) than the 
NLO ones, and rates at LHC energies can be predicted with an accuracy of about 
5%, whereas at Tevatron they are more than 2a above the NLO ones. The NNLO 
predictions for SM Higgs boson production via the dominant gluon fusion process 
have a total (pdf and scale) uncertainty of about 10% at LHC which almost doubles 
at the lower Tevatron energies; they are typically about 20% larger than the ones at 
NLO but the total uncertainty bands overlap. 



1 Introduction 



Parton distributions and their implications have been recently studied within the dynamical 

(radiative) parton model approach up to next-to-next-to-leading order (NNLO) of QCD 

[1]. Here the predicted steep small Bjorken-x behavior of structure functions is mainly due 

to QCD-dynamics at a; ^ 10~^, since the parton distributions at ^ 1 GeV^ are QCD 

radiatively generated from valencelike positive definite input distributions at an optimally 

determined low input scale Qo = /^^ < 1 GeV^. ('Valencelike' refers to a/ > for all 

input distributions x/(x,/i^) oc x°'f{l — x^f , i.e., not only the valence but also the sea 

and gluon input densities vanish at small Such analyses are usually performed within 

the framework of the so-called 'fixed flavor number scheme' (FFNS) where, besides the 

gluon, only the light quark flavors q = u,d,s are considered as genuine, i.e., massless 

partons within the nucleon. This factorization scheme is fully predictive in the heavy 

quark h = c,b,t sector where the heavy quark flavors are produced entirely perturbatively 

as final state quantum fluctuations in the strong field generated by the initial light quarks 

and gluons. Here the full heavy quark mass mc,6,f dependence is taken into account in 

the production cross sections, as required experimentally [2[ |3l IH |5], in particular, in the 

threshold region. However, even for very large values of Q^, ^ ^cb' these FFNS 

predictions up to next-to-leading order (NLO) are in remarkable agreement [6, 7] with 

deep inelastic scattering (DIS) data and, moreover, are perturbatively stable despite the 

common belief that "noncoUinear" logarithms \n{Q'^/m\) have to be resummed for h = c,b, 

and eventually t. This agreement with experiment even at ^ indicates that there 

is little need to resum these supposedly "large logarithms" , which is of course in contrast 

^Alternatively, in the common "standard" approach the input scale is fixed at some arbitrarily chosen 
Qo > 1 GeV^ and the corresponding input distributions are less restricted. For example, the observed 
steep small- a; behavior (a/ < 0) of structure functions and consequently of the gluon and sea distributions 
has to be fitted. Furthermore the associated uncertainties encountered in the determination of the parton 
distributions turn out to be larger, particularly in the small- cc region, than in the more restricted dynamical 
radiative approach where, moreover, the evolution distance (starting at Qq < 1 GeV^) is sizably larger 
(see, e.g., [1 and references therein.) 
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to the genuine coUinear logarithms appearing in hght (massless) quark and gluon hard 
scattering processes. It should be mentioned that, so far, the heavy NNLO 0{a^) 3-loop 
corrections to F2 ^ have been calculated only asymptotically for ^ ml [H [9l [101 E] ■ 

In many situations, calculations within this factorization scheme become unduly compli- 
cated (for a recent discussion, see [12] )• Thus it is advantageous to consider the so-called 
"variable flavor number scheme" (VFNS) despite the somewhat questionable resumma- 
tions of heavy quark mass effects using massless evolution equations, starting at unphysical 
"thresholds" = m^. Here the heavy quarks (c, b, t) are considered to be massless partons 
within the nucleon as well, with their distributions h{x, Q^) = h{x, Q"^) being generated, 
up to NLO, from the boundary conditions h{x,mf^) = h{x,m1) = 0, and at NNLO from 
h{x,mD = h{x,mD = 0{al) as will be exphcitly given in the next Section. Thus this 
factorization scheme is characterized by increasing the number of flavors nf of massless par- 
tons by one unit at = ml starting from n/ = 3 at = Hence the nj > 3 "heavy" 
quark distributions are perturbatively uniquely generated from the Uf — l ones via the mass- 
less renormalization group evolutions (see, e.g. [131 E]; a comparative qualitative and 
quantitative discussion of this zero-mass VFNS and the FFNS has been recently presented 
in [I2j). Eventually one nevertheless has to assume that these massless "heavy" quark 
distributions are relevant asymptotically and that they correctly describe the asymptotic 
behavior of DIS structure functions for scales ^ m^. However, for most experimentally 
accessible values of Q^, in particular around the threshold region of heavy quark (hh) pro- 
duction, effects due to finite heavy quark masses mh can not be neglected. One therefore 
needs an improvement of this zero-mass VFNS where heavy quark mass-dependent correc- 
tions are maintained in the hard cross sections. Such improvements are generally referred to 
as the general-mass VFNS and there exist various different model-dependent ways of imple- 
menting the required mh dependence [I5l|16l[17l[l8l|l3[20linil23[23l|2a|2a 

•^Notice that it is rather superfluous to argue about the best choice of a factorization scheme since the 
scheme choice remains merely a theoretical convention as long as there are no observable signatures which 
allow to uniquely distinguish between the FFNS and any version of a general-mass VFNS (except the 
strictly massless VFNS which has been known to be experimentally inadequate for a very long time.) 
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These factorization schemes interpolate between the zero-mass VFNS (assumed to be cor- 
rect asymptotically) and the (experimentally required) FFNS used for our previous analysis 

In order to avoid any such model ambiguities we shall generate in the next Section the 
"heavy" zero-mass VFNS distributions using our unique NNLO dynamical FFNS distribu- 
tions [1] as input at = m^. This will considerably ease the otherwise unduly complicated 
calculations in the FFNS of gauge- and Higgs-boson production and heavy quark produc- 
tion at collider energies, or the calculation of weak charged-current (anti)neutrino-nucleon 
cross sections at ultrahigh neutrino energies, for example. It has been recently shown 
[12] that for situations where the invariant mass of the produced system {cW, tW, tb, 
Higgs-bosons, etc.) exceeds by far the mass of the participating heavy flavor, the VFNS 
predictions deviate rather little from the FFNS ones, typically by about 10% which is 
within the margins of renormalization and factorization scale uncertainties, and ambi- 
guities related to presently available parton distributions. Let us consider, for example, 
hadronic production. The relevant heavy quark contributions at LO have to be cal- 
culated via gs{d) — cW^, gu — hW^ in the (fully massive) FFNS as compared to the 
much simpler quark fusion subprocesses cs{d) W~^, bu in the VFNS, etc. Here 

nonrelativistic contributions from the threshold region in the FFNS are suppressed due 
to ^/hh/'fnc,b — Mw/rric^b ^ 1. Similarly, hadronic single top production via VT-gluon 
fusion [29] requires in the FFNS the calculation of the subprocess ug dtb at LO and of 
ug —>■ dtbg, etc., at NLO; in the VFNS one needs merely ub — > dt at LO and ub —>■ dtg, 
etc., at NLO, using massless initial-state partons. Again, ^/su^/mh — mt/rrih ^ 1 and thus 
the FFNS and VFNS results are not too different. A similar agreement is obtained for 
hadronic (heavy) Higgs boson production where the LO FFNS subprocess gg bbH has 
to be compared with the bb fusion subprocess (for massless initial-state partons) in the 
VFNS starting with bb ^ H at LO. (H = H^^; h^, H^, A° denote the Standard Model 
(SM) Higgs boson or a light scalar a heavy scalar and a pseudoscalar of super- 
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symmetric theories with Mh ~ 100 GeV.) Again, ^fsth/'fnh = (Sm^ + MH)/'mb ^ 1 in the 
FFNS which indicates that the simpler LO, NLO, and NNLO VFNS bb fusion processes 
do provide rehable predictions. (Notice that these situations are very different from DIS 
heavy quark hh production via 'y*g — hh, etc., where \/Sth/mh = 2 is not sufficiently large 
to exclude significant contributions from the threshold region and therefore the VFNS 
predictions deviate sizably from the FFNS ones |12j). 

Within the present intrinsic theoretical uncertainties we can therefore rely on our 
uniquely generated NNLO VFNS parton distribution functions (pdfs) where, moreover, 
the required NNLO cross sections for massless initial-state partons are, in contrast to the 
fully massive FFNS, available in the literature for a variety of important production pro- 
cesses. The perturbative stability of the NNLO predictions, when compared with the ones 
based on our dynamical NLO VFNS pdfs [12], will be furthermore studied in the next 
Section for the hadronic production of and bosons, as well as of the SM Higgs 
boson at the Tevatron and at LHC. Our conclusions are summarized in Sect. 3. Finally, 
the Mellin n-moments of the renormalized heavy-quark flavor operator matrix elements 
relevant for the generation of the VFNS pdfs at NNLO are summarized in the Appendix. 



2 Heavy flavor parton distributions and their impli- 
cations at high energy colliders 



As common, the flavor transitions Uf Uf + 1 are made when the factorization scale 
equals the (pole) mass of the heavy quarks, = ml, and the pdfs for Uf + 1 flavors are 
defined from the light flavor pdfs and the massive operator matrix elements for nf light 
flavors. In Mellin ra-moment space, the heavy quark pdfs can then be expressed in terms 
of the original light ones at NNLO as 
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{h — h)nj:+i{n,mD = and the remaining matching conditions for the hght pdfs and the 
gluon distribution read 

(g ± {n, ml) = {q ± g)„^. {n, ml) + A^^f^ (n) (g ± g)„^ {n, ml) (2) 

gnf+i (n, ml) = Qn^ {n,ml) + al 
with the moments of the flavor singlet quark distribution being given by 

Sn,(n,g2) = / dx V [qkix,Q'') + qk{x,Q')] (4) 

where qi = u, q2 = d, etc. The coefficients A^'^\x) of the operator matrix elements have 
been originally calculated in [TB] and their Mellin moments A^'^\n) have been analyzed 
and given in [HI HU |30]. Due to our choice = ml for the thresholds, only the scale- 
independent parts of the expressions for A^'^(n) are needed which, for completeness, will 
be summarized in the Appendix. The strong coupling = as{Q'^)/4:iT is matched at the 
various thresholds = ml in the standard way as recapitulated in pj with q;s(M|) = 
0.1124 as obtained in our dynamical scenario using = 1.3 GeV, = 4.2 GeV, and 
mt = 175 GeV. Our choice for the input of the 'heavy' VFNS distributions in (1) are the 
unique NNLO dynamical FFNS distributions [T] at = m^, as obtained from the NNLO 
evolution of our valencelike input distributions at Q"^ = ji'^ = 0.55 GeV^ (see Table I of 
[1]. The resulting VFNS predictions at scales ^ ml should become insensitive to this 
input selection [13], since asymptotically the VFNS pdfs are dominated by their radiative 
evolution rather than by the specific input at = ml, i.e., because of the long evolution 
distance input differences get evolved away at ^ m| where the universal perturbative 
QCD splittings dominate. 

For illustration we show in Fig. 1 our NNLO charm and bottom distributions together 
with the fully convoluted and F2 structure functions which are also compared with the 
NLO ones. In general the NNLO results for Fg'^ fall below the NLO ones (dash-dotted 
curves) in the small-x region. Here at NNLO the 0{a'^) convolutions of the fermionic and 
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gluonic coefficient functions with h and the gluon distribution, respectively, become more 
important than at NLO since the 'heavy' quark distributions xc and xb by themselves 
(short-dashed curves) are sizably different from and respectively. At NLO the 
0{as) the quark and gluon convolution contributions almost cancel |12j and thus xc and 
xb almost coincide with the appropriate NLO structure functions in Fig. 1. As is obvious 
from Fig. 1, however, such differences between NNLO and NLO results lie always within 
the la — 2a uncertainty bands in the relevant large region, >, 10^ — 10'^ GeV^, and 
can therefore be hardly delineated experimentally. 

The shape of the gluon distribution at two typical fixed values of x, relevant for Higgs 
boson production at LHC, is illustrated in Fig. 2. At small to medium values of the 
NNLO gluon falls always below the NLO one and in both orders the gluon remains positive 
at small in the very small-x region. This dampening of the NNLO gluon is a typical 
NNLO effect being mainly caused [1] by the gluonic 3-loop splitting function Pgg which is 
negative and more singular (~ — Mn ^) in the small-x region [31J than the NLO (and LO) 
ones. At large values of the NNLO and NLO gluon distributions become practically 
indistinguishable . 



2.1 Weak gauge boson production 



As a next test of our VFNS distributions we turn to the hadronic and production. 
The inclusive differential cross section is usually written as |32] 

da^ 

— = rav{Q',M^)Wv{T,Q'), r = / s (5) 

where V is one of the gauge bosons of the Standard Model (7, or W^) which subse- 
quently decays into a lepton pair {I1I2) with invariant mass M^^i^, i.e. = M^^^^, and 
ay is the pointlike cross section, e.g., a^ = Ana"^ /9Q^, etc. [32j. The hadronic Drell-Yan 
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structure function is represented by 




r 

t/xi ^2 



Xi,X2 



T 



(6) 



with -P-Dj^ denoting the usual combination of pdfs of flavor type % and j which depend on 
the factorization scale /ii?. The QCD correction term is expanded in a power series of as 
(or agj^TK or a^/Tr) as follows 



with ^Xj and the NNLO 2-loop A^^ being given in [321 133] , and the choice for the renor- 
malization scale [i^ = fip is dictated by all presently available pdfs which have been 
determined and evolved according to fiR = fip- The scale uncertainties of our predictions 
are defined by taking Mvy/2 < /xf < 2Mw, using Mw = 80.4 GeV (and similarly for 
production, using Mz = 91.2 GeV). Furthermore, it should be noted that only the initial 
u,d,s,c quark flavors and the gluon contribute sizably via the various fusion subprocesses 
in (6) to the production rates of gauge bosons, whereas all subprocesses involving the 
6-flavor distribution, e.g., ub — > W~^, cb — > W~ , etc., are negligibly small [T2] . 

Our NNLO predictions for a{pp —>■ W^X) and a{pp — * Z^X) are compared with our 
NLO ones [12] in Fig. 3 where, for comparison, we also show the predictions of Alekhin 
[HI IM] . The vector boson production rates at NNLO are typically slightly larger (by more 
than la than at NLO with aK = NNLO/NLO factor of K^^'+w- = i.o4 and K^' = 1.06 
at Tevatron energies {^/s = 1.96 TeV, cf. Table 1), to be compared with the predictions 
of Alekhin [HJ |31] K^^~^^ ~ = 1.03. This confirms again the fast perturbative 
convergence at NNLO since the NLO/LO i^-factor [12] is 1.3 for + W~ production 
at y/s = 1.96 TeV. Our predicted NNLO cross sections at ^/s = 1.96 TeV (cf. Table 1), 
a{pp —>■ + W~ + X) = 25.2 nb and a{pp ^ Z^ + X) = 7.5 nb, are similar to the 
ones of MSTW t35j, 25.4 nb and 7.4 nb, respectively, but smaller than the ones obtained 
by Alekhin pjj, 25.8 nb and 7.8 nb, respectively. For the latter cases the branching ratios 
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B{W iu) = 0.108 and B{Z £+£-) = 0.034 have been used. It is obvious from Fig. 3 
that most of these results are within the present experimental Icr uncertainty. The scale 
uncertainties of our NNLO predictions in Fig. 3 at ^/s = 1.96 TeV amount to less than 
0.5% (where fip = My/2 gives rise to the upper limits and fip = 2My to the lower limits, 
with V = W^, which is four times less than at NLO [12J. 

In Table 2 we present our NNLO predictions for and Z° production at LHC energies. 
For comparison we also display our previous NLO results [12j. Here the scale uncertainties 
amount to less than 1.7%, i.e., are about half as large than the stated pdf uncertainties 
and than the scale uncertainties at NLO [I2]. For example, the full NNLO expectations at 
a/s = 14 TeV are 

a{pp^W+ + W- +X) = 190.2 ±5.6pdf ^};^|scaie n6 (8) 
+ = 55.7±1.5pdf t'oiUie nb . (9) 

Here the scale choice = 2My gives rise to the upper limits, and fxp = My /2 to the lower 
limits of our predicted cross sections. These results are about 5% smaller than the ones 
of MSTW [25], whereas Alekhin [M] obtained 195.2 nh and 57.7 nh for W+ + W~ and Z° 
production, respectively, with similar pdf uncertainties as in (8) and (9). /^From Table 2 
it becomes obvious that the vector boson production rates somewhat increase at NNLO 
as compared to the NLO expectations, but such differences are well within present pdf 
and scale uncertainties. Moreover, the smallness of such differences {K ~ 1.02) indicates 
the reliability of perturbative predictions already at NLO. For comparison we note that 
within the FFNS (where the heavy c, 6, t quark flavors do not form massless partons of the 
nucleon) the + W~ production rate has been estimated [12] to be about 192.7 nh at 
NLO with a total (pdf as well as scale) uncertainty of about 5%. In general the NLO- 
VFNS prediction of 186.5 nh in Table 2 falls somewhat below that estimate but remains well 
within its total uncertainty of about 6% [12]. Due to the reduced scale ambiguity at NNLO 
and due to the slightly different NNLO estimates obtained by other groups as discussed 
above, we conclude that the rates for gauge boson production at LHC energies can be 

8 



rather confidently predicted with an accuracy of about 5% irrespective of the factorization 
scheme. 

2.2 Higgs boson production 

As a final application of our NNLO VFNS pdfs we consider the hadronic production of 
the SM Higgs boson. Similar as for gauge boson production in Sec. 2.1, the total inclusive 
cross section for Higgs boson production is usually written as [32|, |33l HOI HI] 

= / dxi dx2 PDfj (xi, X2, f/p) aij^n (s = xi X2 s, jjp, /i^) (10) 
with the partonic cross sections for ij HX being written as 

2 

aij^His) = aoA,j{s) = ^ < (/i^) A[f{s, /i^) (11) 

n=0 

where the obvious fip and dependencies have been suppressed. The dominant Higgs 
production proceeds via gluon-gluon fusion where uo = /{57Qnv^) for the initial LO 
process gg H, with the Higgs vev v = {V2Gf)-^/^ ^ 246 GeV, and the NLO QCD 
corrections A^-j'* are given in [321 113] and the NNLO A\^^ ones in [33l HQ] . The factoriza- 
tion scale is ususally chosen to be fip = Mh, and the scale uncertainty is illustrated by 
taking < fip ^ 2Mh- The much smaller contribution stemming from bottom quark 

annihilation starts at LO with hh H where (Tq = 7rA^/(12M^) with A;, = \/2m}j/v in 
the SM, and the NLO and NNLO QCD corrections Ag^ and Af^ can be found in g!]. 
Here it has been argued [HI HSj HH] that the optimal choice of the factorization scale 
is lip — Mh/4: where the differences between the VFNS and the FFNS are significantly 
reduced and the LO, NLO and NNLO results become rather similar [H], which implies 
a more stable perturbative behavior. The scale uncertainty is again probed by taking 
^{Mh/4:) < fip < 2{Mh/A). In both cases ^ir = ^ip as dictated by all presently available 
pdfs. 
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In Fig. 4 we show, as a function of the Higgs mass, our NNLO (thick sohd curve) and 
NLO (thick dash-dotted curve) predictions for LHC for Higgs boson production via the 
dominant gluon-gluon fusion subprocess which starts, at LO, with gg H . The shaded 
regions around these central predictions are due to the ila pdf uncertainties. Reducing 
the scale to \Mh one arrives at the thin upper curves at each order, whereas the scale 
choice /iF = 2M/f results in the respective lower curves, where the appropriate ±lcr pdf 
ambiguities have also been included for each choice of scale. These ambiguities for each 
scale choice /iF = /^ij = ^Mh, Mh, 2Mh are more explicitly illustrated in Table 3. Despite 
the fact that the NLO and NNLO total uncertainty bands overlap in Fig. 4, the predicted 
NNLO production rates are typically about 20% larger than at NLO. The insensitivity of 
these predictions with respect to the appropriate choice of the pdfs is illustrated by the 
dashed curve which has been obtained by using NNLO matrix elements and (inconsistently) 
NLO pdfs. Here for the dominant gluon fusion process such an inconsistent choice of the 
pdfs appears to be immaterial and the production rates depend dominantly on the NNLO 
QCD dynamics. Our central predictions in Fig. 4 are comparable with the ones presented 
in [17], but are about 10% smaller than the ones in [35]. For completeness we also show 
in Fig. 5 our NNLO expectations for Higgs boson production at the Tevatron, ^/s = 1.96 
TeV. Note that here the total uncertainty bands almost double at NNLO and NLO as 
compared to the ones at LHC in Fig. 4. 

In Fig. 6 we finally show the subdominant contribution to Higgs boson production at 
LHC due to bottom-quark fusion which starts with bb ^ H at LO. Here, in contrast 
to the by far dominant gluon fusion process in Fig. 4, the NNLO and NLO predictions, 
together with their ±la pdf uncertainties, almost coincide with the NNLO results falling 
very slightly below the NLO ones. Here, however, the correct choice of the NNLO pdfs 
turns out to be important, since choosing (incorrectly) NLO pdfs [12] for a NNLO analysis 
results in too small a production rate as shown by the dashed curve. At NNLO the scale 
dependence is here, again in contrast to the by far dominant gluon-gluon fusion process. 
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very marginal: using fip = fJ'R = '^{Mh/4:) instead of Mh/4: leaves the results in Fig. 6 
practically unchanged, whereas the choice [ip = iir = \{Mh/^) increases the results by at 
most 5%. 

It should be again emphasized that here, as in the previous case of gauge boson produc- 
tion, the simpler VFNS yields sufficiently reliable predictions for Higgs boson production 
despite the fact that a fully massive FFNS analysis cannot be performed at NNLO at 
present (due to the absence of NNLO, and in many cases even NLO, matrix elements with 
ruh 7^ 0). This is due to the fact that \/~Sth/^b = (Sm^ + MH)/mb ^ 1, i.e., nonrelativistic 
contributions from the threshold region in the FFNS are suppressed, and thus the FFNS 
and VFNS predictions should not differ too much [12], as has been discussed in more de- 
tail in the Introduction. Indeed it has been noted [12] that the FFNS and VFNS results 
at NLO are compatible [Ml SHI ED], and that the VFNS rates exceed the corresponding 
FFNS Higgs boson production rates by about 10-20%, depending on the choice of the scale 



3 Summary and Conclusions 



Based on our recent NNLO dynamical parton distributions as obtained in the FFNS pQ, 
we generated radiatively VFNS parton distributions at NNLO where also the heavy quark 
flavors (c, b, t) become massless partons within the nucleon. The latter pdfs in the 'variable 
flavor number' factorization scheme considerably ease the otherwise unduly complicated 
calculations in the FFNS where for the time being fully massive NNLO analyses are not 
possible (and in many cases even not at NLO), such as the calculation of gauge- and 
Higgs-boson production and heavy quark production at collider energies. It has been 
shown [12] that for situations where the invariant mass of the produced system exceeds 
by far the mass of the participating heavy flavor in the FFNS, the VFNS predictions 
deviate rather little from the FFNS ones, typically by about 10% which is within the 
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margins of renormalization and factorization scale uncertainties and ambiguities related 
to presently available parton distributions. As an application of our NNLO VFNS pdfs 
we studied the perturbative stability of the predictions for gauge (W"^, Z^) and SM Higgs 
boson production at collider energies by comparing them with the appropriate NLO results, 
taking into account pdf uncertainties as well as scale dependencies. The NNLO predictions 
for gauge boson production are typically slightly larger (by more than lo") than the NLO 
ones, cf. Table 1 and 2. Due to the reduced scale ambiguity at NNLO and due to the 
slightly different NNLO estimates of other groups we conclude that the rates for gauge 
boson production at LHC energies can be rather confidently predicted with an accuracy 
of about 5%. At the Tevatron {^/s = 1.96 TeV) the NNLO predictions are more than a 
2cr pdf uncertainty above the NLO ones, but most of these results are within the present 
experimental la uncertainty. 

The NNLO predictions for the production of the SM Higgs boson via the dominant 
gluon fusion (in contrast to the subdominant bottom-quark fusion) process are, at collider 
energies, typically about 20% larger than at NLO, but their respective total (pdf and 
scale) uncertainty bands overlap. Higgs boson production at LHC (y^ = 14 TeV) can be 
predicted with an accuracy of about 10% at NNLO (with the total uncertainty being almost 
twice as large at NLO), whereas the uncertainty almost doubles at Tevatron {^/s = 1.96 
TeV). 

A FORTRAN code (grid) containing our NNLO- VFNS pdfs (including their uncer- 



tainties) can be obtained on request or directly online from |http:/ /doom.physik.uni- 



dortmund . de/p df server 
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Note added. 

While completing this manuscript, an investigation along similar lines appeared [S. Alekhin, 



J. Bliimlein, S. Klein, and S. Moch, arXiv: 0908.2766]. Several results are similar to ours. 
However, the gauge boson production rates are about 4% larger at the Tevatron and about 
10% larger at LHC than our NNLO predictions. Similarly, the NNLO predictions for Higgs 
boson production at LHC {y/s = 14 TeV) are 5-8% larger for ^ 150 GeV than ours, 
but agree with us for larger Higgs masses, whereas at the Tevatron their expected rates are 
12-30% smaller than our ones for Mh =100-200 GeV. The comparison of these production 
rates refers always to the central results, disregarding all pdf and scale uncertainties. 
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Appendix 



Some analytic form of the Mellin n-moments of the operator matrix elements including 
the heavy quark flavors, which have been originally calculated in Bjorken-x space [16] and 
which are needed in (l)-(3), have been already implicitly used in some NNLO evolution 
programs (see, e.g. [51]). Here we summarize the relevant analytic expressions which can be 
directly continued to complex values of n as required for the Mellin inversions to Bjorken-x 
space. 

The Bjorken-x expressions are given in Appendix B of [16j and we follow the notation 
used there [cf. also (l)-(3)]. Due to our choice = ml for the flavor transition thresholds, 
only the scale-independent parts of these expressions contribute. Their moments are as 
follows: 

CpTf ^'^ ^""^ ~ (n-l)n2(n + l)2(n + 2) '^"^ ' 27 n - 1 n ^ 

24 12 56 28 24 1960 1 



n + 1 (n + l)2 {n + If {n + lf 27 n + 2 
48 1 64 1 

9" (n + 2)2 + y {n + 2f 



448 1 64 1 ,^ ^, 



where = |, T/ = | und Sk{n) = YJj=i jk using 

5i(n)=^(n + l)+7i., 5,,(n) = |^;^#^'-^)(n + l) + C(fc'), k' > 2, (A.2) 

with = S'+^HnT{z)/dz'+^ and 'Je = 0.5772156649, C(2) = and C(3) = 

1.2020569032 for the analytic contiunation to complex n. Since the moment of the rather 
complicated coefficient A^'^ appearing in (1), and as given in (B.3) of ^16j, cannot be 
straightforwardly expressed in terms of analytic functions of n [8], we have employed for 
practical purposes the ra-moment of the sufficiently accurate x-parametrization suggested 
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in [51]: 

75f2). X ....S'fin) ^ ,S?(n) ,2.77 293.6,^,, 3.333^,,^,, 

Alf{n) = 1.111^^-0.4^^+ + ^^)S,{n) + S,{n)S2{n) 

^ n n n v? n 

,0.4 293.6, ^ , , S-i(n) 24.89 
-( —) S2(n) + 295.822 - 0.006 

^ n n n — 1 

187.8 + 293.6C(3) , 93.68 - 293.6C(2) 6.584^9.336^249.6 ^^^^ 



The remaining coefficients relevant for the hght quark and gluon sector in (2) and (3) 
respectively, can be straightforwardly transformed to n-space : 

Koh K^) = -^Si{n-1) + — S2{n-1)--S3{n-1) + — + 



C^Tf '''^^'^ ~ 27 " 9 ^^^'^ 3 "^'^'^ " 18 " 27 n 9 

268 1 44 1 4 1 4 1 

n + 1 ^ "g" (n + 1)2 " 3 n3 " 3 (n + 1)3 



3 n - 1 
8 1 



1 56 

Sl{n - 1) - Si (n - 1) + S2{n - 1) + - 



+ 



3 n 
4 



3^[5?(n+l)-f 5i(n + l) + 5^(^ + 1) + 



(A.5) 



4CpT, 



/ 



15 



2 20 

n — 1 n 



+ 



+ 
6 



12 



12 



n 



(n + l)2 



+4C^r, 



(n + l)3 (n + l)4 n + 2 

56^, lSi(n + l) 5 139 1 

27^^(^-^)+3^Tr^+ 18+^;r3T 



157 



22 



1 



+ 



1 



27n 
175 1 



2 137 1 

^3n3^^^^^1 9 (n + l)2 ' 3(n+l)3 27 n + 2. 



13 
(A.6) 



with Ca = 3. 
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^Pi5--vx (j^^)^ ^=1.96 TeV 


V 


NNLO 


NLO 


W+ + W" 

z° 


12.6 ±0.1 
25.2 ±0.3 
7.5 ±0.1 


12.1 ±0.1 

24.2 ±0.3 
7.1 ±0.1 



Table 1: NNLO predictions for vector boson production at the Tevatron, with the NLO 
ones being taken from ^2]. The errors refer to the ±lcr uncertainties imphed by our 
dynamical NNLO [T] and NLO pdfs. The scale uncertainties of our NNLO predictions, 
due to ^My < /Ui? < 2Mv, amount to less than 0.5% (i.e., are about half as large as the 
stated pdf uncertainties) which is about four times smaller than at NLO [T2] . 
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^pp-vx ^=10 TeV 


V 


NNLO 


NLO 


W+ 
W+ + W" 


78.7 ±2.0 

55.8 ± 1.4 
134.5 ±3.4 
39.1 ±0.9 


76.7 ± 1.7 
54.7 ± 1.2 
131.6 ±2.9 
38.1 ±0.8 




^pp-vx (^^^^^ 7^=14 TeV 


V 


NNLO 


NLO 


w+ 

W+ + 

z° 


109.8 ±3.2 
80.4 ±2.4 
190.2 ±5.6 
55.7 ± 1.5 


107.5 ±2.9 
79.1 ± 2.1 
186.5 ±4.9 
54.6 ± 1.3 



Table 2: As in Table 1 but for LHC energies. The scale uncertainties of our NNLO 
predictions amount to less than 1.7% of the total predicted rates which is about half as 
large as the stated pdf la uncertainties and the scale uncertainties at NLO [12j. 
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\Mh 


Mh 


2Mh 


100 


67.6 ±3.0 


62.2 ±2.6 


57.3 ±2.2 


150 


33.0 ± 1.2 


30.4 ± 1.0 


28.1 ±0.9 


200 


19.7 ±0.6 


18.3 ±0.5 


16.9 ±0.5 


250 


13.6 ±0.4 


12.6 ±0.4 


11.7 ±0.3 


300 


10.7 ±0.3 


9.9 ±0.3 


9.2 ±0.3 



Table 3: Typical NNLO scale dependencies of the cross sections (in units of pb) for Higgs 
boson producion at a/s = 14 TeV via the dominant gluon-gluon fusion subprocess with 
Mh in GeV units. The errors refer to the la pdf uncertainties. The maximal upper limits 
at Hf = IJ'R = \Mh agree with the thin solid curve at NNLO in Fig. 4, whereas the lower 
curve in Fig. 4 corresponds to the minimal lower limits at /Xi? = /xr = 2Mh- 
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Figure 1: The predicted x-dependencies of the charm and bottom quark structure func- 
tions ^F2{x,Q^) and respectively, in the zero-mass VFNS, together with their 
±l(j uncertainties, at some typical fixed values of Q^. The NNLO charm and bottom 
distributions, xc{x,Q'^) and xb{x,Q^), are shown by the short-dashed curves. 
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Figure 2: The NNLO and NLO gluon distributions together with their ±la uncertainty 
bands at two representative fixed values of x. 
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Figure 3: Predictions for the total + W~ and production rates at pp colliders with 
the data taken from [SUl EH EHl EH]- Our NLO VFNS predictions are taken from [12], 
and the NLO and NNLO ones of Alekhin from [TU The adopted momentum scale is 
fJ^F = fJ'R = My for V = W^, Z^. The scale uncertainties of our NNLO predictions, due 
to \Mv < /^f < IMy, amount to less than 0.5% at y/s = 1.96 TeV, i.e., is four times less 
than at NLO ^2]. The shaded band around our NNLO and NLO predictions are due to the 
±la uncertainty implied by our dynamical NNLO [T] and NLO [7] parton distributions. 
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Figure 4: Predictions for SM Higgs boson production at LHC {pp — » HX) via the dominant 
gluon-gluon fusion process, which starts at LO with gg ^ H via a top-quark loop. The 
shaded bands around the central NNLO and NLO predictions are due to the ilex pdf 
uncertainties, all referring to a scale choice fip = fJ^R = Mh- The thin solid and dash- 
dotted curves above these NNLO and NLO bands refer to a scale = fJ'R = with 
±lcr pdf uncertainties included, and similarly the lower curves refer to fip = fJ'R = ^Mh 
(for more details cf. Table 3). The dashed NNLO curve is obtained by using NNLO matrix 
elements and (inconsistently) NLO pdfs [12] with iip = 1^r = ^h- 
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Mh (GeV) 

Figure 5: As in Fig. 4 but for the Tevatron [pp —>■ HX). 
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Figure 6: Predictions for SM Higgs boson production at LHC via the small subdominant 
bottom-quark fusion process which starts with bb H at LO. The shaded bands corre- 
spond to the ±lcr pdf uncertainties of the NNLO and NLO central predictions, all referring 
to a scale choice fip = fJ^R = Mh/4:. The dashed NNLO curve is obtained by using NNLO 
matrix elements and (inconsistently) NLO pdfs [12]. 
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